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Abstract 

We show that a recently proposed extension of the minimal supersym- 
metric standard model can provide a scenario where both the cold and the 
hot dark matter of the universe owe their origin to a single scale connected 
with the breakdown of global B — L symmetry. The susy partner of the 
majoron and the light Majorana neutrinos are the cold and hot dark matter 
candidates respectively in this model and their desired relative abundances 
emerge when the scale of B — L symmetry breaking is in the TeV range. 
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The idea that the large scale structures seen today evolved from very 
small primordial density inhomogeneities has been strengthened by the recent 
detection of large scale anisotropics in the cosmic microwave background ^ . 
Nevertheless, one of the necessary ingredients for the structure formation, 
namely the nature of dark matter, remains unknown. 

The most satisfactory model for structure formation is perhaps the cold 
dark matter (CDM) theory 0] where the Universe is assumed to be spatially 
fiat {Q = 1) and with ~ 0.9 of the mass density formed by CDM particles, 
cold in the sense that they decouple from the expanding thermal bath at 
temperatures much smaller than their mass. CDM can successfully explains 
galaxy-galaxy and cluster-cluster correlation functions on scales of order of 1- 
5 Mpc. However, it now appears to be inconsistent with large scale structure 
data in the automatic plate machine (APM) galaxy survey f^, which suggest 
more power on large scales than the standard CDM predictions. On small 
scales, the observed pairwise velocity dispersion for galaxies appears to be 
smaller than those predicted by CDM IQ. 

One alternative is the hot dark matter model (HDM). HDM is taken 
to be a light neutrino, which decouples from the thermal bath when still 
relativistic, with ttIj, = (92 Q^, K^) eV where H = 100 h Km/sec Mpc is 
the Hubble parameter. In the pure HDM picture, the processed fluctuation 
is characterized by the distance a neutrino travels over the history of the 
Universe, Aj, ^ 40 (30 eV/m^) Mpc. The problem with this however is that 
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Aj, is too large with respect to the scale which is just becoming nonlinear 
today, ~ 5/i^^ Mpc. If galaxy formation occurs early enough to be consistent 
with high-redshift galaxies and quasars, then structures on 5 Mpc will 
overdevelop. 

The hope is that cold + hot dark matter (C + HDM) will combine the 
success of both models. Indeed (C + HDM) models with VLcdm — 0.6, 
Vty ~ 0.3, flbaryon — 0.1 and a Hubble constant h ~ 0.5 provide the best fit 
for microwave anisotropy data, large scale structure surveys, and measures 
of the bulk flow with a few hundred megaparsecs [Q. 

Even if C + HDM is appealing for the large scale structure phenomenol- 
ogy, it might seem rather unpalatable from the point of view of particle 
physics for the following reasons. It is well-known that low-energy supersym- 
metric theories provide an elegant solution to the hierarchy problem [^. The 
lightest supersymmetric particle (LSP) is the most attractive candidate for 
CDM and is made stable by imposing a discrete symmetry, called i?-parity. 
Nevertheless, in the minimal supersymmetric standard model (MSSM), neu- 
trinos are not massive and therefore there are no candidates for HDM. If 
one considers, for instance, minimal extension of the MSSM by including a 
right-handed neutrino superfield u'^ and a singlet field S with lepton number 
L{S) = +2 [|^, the U{1)b~l and i?-parity symmetries get spontaneously bro- 
ken, neutrinos acquire a mass, but the LSP is no longer stable and cannot play 
the role of CDM. So, when trying to build up a supersymmetric model with 
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both CDM and HDM components, one has to face the problem of making 
neutrinos massive while at the same time, keeping the LSP stable. Moreover, 
even if one is successful in doing so, the relative abundances of light mas- 
sive neutrinos and LSP's are expected to be set by two uncorrelated scales, 
namely the U{1)b-l breaking scale, Vbl, and the supersymmetric breaking 
scale, Ms, and a sort of fine-tuning on Vbl and Ms must be done to satisfy 
the requirements of C + HDM model. 

In this Letter, we give an example of an extension of the MSSM, where it 
is possible to implement C + HDM model for large scale structure formation 
without arbitraily separating the scales Vbl and Ms- The cold dark matter 
in this model is the lightest supersymmetric particle of the theory ( not the 
usual neutralino ) with mass in the range of 10 to 50 GeV while the hot dark 
matter is an eV Majorana neutrino. What we find quite interesting is that 
the relative abundances of the HDM and CDM components are set by the 
same scale Vbl — 10^ GeV. The new scale being so low brings the hope that 
the effects of this low scale may be testable in near future. 

The model is a supersymmetric version of the original singlet Majoron 
model and was recently proposed in ref. 0. For Vbl > Ms (where Ms 
is the scale of supersymmetry (SUSY) breaking), the susy partners of the 
majoron (the smajoron and majorino) acquire mass of Ms . Since they are 
very weakly coupled to light particles of the theory such as the quarks and the 
leptons, for their cosmological relic abundance not to disturb nucleosynthesis. 
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the scale oi B — L symmetry breaking Vbl must have an upper bound of order 
of a few TeV§. In this letter, we show that if the majorino is the lightest 
superparticle of the model (i.e. lighter than the familiar neutralino), then 
due to i?-parity conservation the majorino is stable and for Vbl — TeV, has 
the right relic abundance to be the cold dark matter. The neutrinos in this 
model are of course massive with either the or or both having mass in 
the eV range and therefore becoming the obvious hot dark matter candidate. 
As far as the Ur is concerned, cosmological constraints demand that in our 
model its mass either be in the eV range ( thereby sharing the dark-matter- 
dom with other neutrinos ) or in the MeV range in which case it decays to 
i/e,/i + Majoron fast enough so as not to effect structure formation in the 
universe. 

Let us now proceed to discuss the details of the model. The superpotential 
can be written as the sum of two terms 

W = Wo + Wi, (1) 

where Wq is the usual MSSM piece and (generation indices are suppressed) 

Wi = KLH2V^ + fv^v'^Si + A (S1S2 - M^) Z. (2) 

Here M is an explicit mass scale of order of Vbl ; is a right-handed neutrino 
superfield which is a singlet under SU{2)l®U{1)y and carries lepton number 
L = —1; Si, S2 and Z are SU(2)l ® U{1)y singlet superfields as well and 
carry lepton numbers L = 2, —2 and 0, respectively. 
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It easy to work out from eqs. (1) and (2) the full scalar potential p. In 
particular, we are interested in the pieces containing the fields Si, S2 and 
Z (we are assuming h,^ <^ A, / and all the couplings real for simplicity) 



V = 4f\9''f\Sif + X^SiS2-M''' + X'\Si\'\Z 

+ \fi)'i)^ + XS2Zf + ml \Sif + ml l^al^ + m| \Zf + mj. 

+ (AfD'u'Si + h.c.) + (AxSAZ + h.c.) - [BxM^Z + h.c.) . (3) 

In eq. (3) {mi,m2,mz,m^c, Aj, Ax, Bx) are all parameters which set the 
scale of the soft breaking of supersymmetry and we assume them to be of 
the same order of magnitude, Ms — 10^ GeV. 

One would naively expect that for temperatures below Ms, the symmetry 
breaking effects would be negligible in analogy with what happens in spon- 
taneously broken gauge symmetries such as electroweak. However, the above 
symmetry breaking is not spontaneous but explicit and there is no reason to 
expect the restoration of the symmetry above Ms- Rather, the relevant ques- 
tion is above which temperature, the effect of the soft symmetry breaking 
terms falls out of equilibrium. One could estimate easily this temperature by 
comparing the rate for the symmetry breaking effects Fg ~ M|/T with the 
expansion rate of the Universe, H ~ g^,(TY^'^T'^ / Mpi, where g^.{T) counts 
the effective number of relativistic degrees of freedom and Mpi ^ 1.2 x 10^^ 
GeV is the Planck mass. One finds that the symmetry breaking effects are 
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out of equilibrium for 

TZTs^"^^^^ ^IQ'GeV. (4) 

Since we are interested in values of Vbl not far from Ms (see later) and 
certainly in the range Ms^Vbl^Ts, we cannot neglect in the scalar potential 
(3) the soft supersymmetric breaking terms. 

As the temperature falls below the value T ^Vbl, the U{1) b-l symmetry 
is spontaneously broken by the following vacuum expectation values (VEV's) 

{S^) = V^, (^2) = ^2, V^^V^^Vbl, (5) 
whereas it is easy to show that Z acquires a VEV 

= x^v,^ + vi) ■ 

The crucial point here is that, whenever Vbl^Ms, the right-handed sneutrino 
D'^ docs not acquire a VEV so that the discrete i?-parity symmetry, which is 
proportional to (—1)^, is preserved. As a consequence, the LSP in our model 
remains stable and can provide a suitable candidate for CDM. 

After the spontaneous breaking of the U{1)b-l symmetry, a Nambu- 
Goldstone boson, the majoron, will appear 

J ^ V,{lmSi)-V2{lmS2) 

whereas its fermionic superpartner, the majorino 

V'J = ; (8) 
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and its real superpartener, the smajoron 



V,{ReS,)-V2iReS2) 

(^J = ; (9) 



will acquire a mass proportional to Ms- In particular, the majorino ipJ 
receives two different mass contributions: one, at the tree level, of order of 
X{Z), and the second from one- loop diagrams involving u'^ and u'^ of order 
of {fMs/lQir'^). Without any fine-tuning of the parameters, we can have 
m^j ^ (10 — 50) GeV and consider the majorino lighter than any other 
neutralino. As a consequence, the majorino will be the LSP and a suitable 
candidate for CDM. The smajoron being heavier can decay through various 
channels as discussed in ref. 0. We have found a new decay channel aj 
J J [0, which provides a somewhat faster decay mode. We find that, this 
decay amplitude goes like ~ (Mg/V^j^) and requiring the aj lifetime to be 
~ 10^^ sec. in order not to disturb neucleosynthesis discussion implies that 
Vbl < 10^ GeV. Our discussion of majorino as the CDM needs Vbl — TeV, 
which is consistent with the above bound. 

We now proceed to the estimation of the relic abundance of majorinos. We 
first want to point out that the couplings of the majorino, the smajoron and 
the majoron among themselves and with the other fields must be calculated 
in presence of the soft supersymmetry breaking terms and not in the exact 
limit of supersymmetry , since we are considering temperatures below T5. 
For T < Vbl all the heavy fields, such as u'^, decay into lighter particles 
so that the majorino, as well as the smajoron and the majoron, can only 



interact among each other and with the "standard model" particles, such as 
leptons and quarks, through the coupling h^LHiu'^. Therefore, the coupling 
of ip'jS and J' s to the "standard model" particles is suppressed by powers of 
e ~ (/^/i^/167r^) and they decouple from the "standard model" thermal bath 
at a temperature [§] 



To get a feeling of the numbers, T^/j ~ 10^ GeV for Vbl — 10^ GeV. There- 
fore, the majorino is expected to decouple from the "standard model" thermal 
bath when still relativistic. Nevertheless, its number density for T<T^d does 
not decrease only due to the expansion of the Universe. Indeed, the key point 
here is that, even after T^d, the number density of majorinos continue to fol- 
low its equilibrium value due to the fact that they keep into equilibrium with 
the thermal bath formed by majorons via interactions mediated by heavy 
particles with mass ~ Vbl such as KeZ or the two fermionic combinations 
of ipi, and ipz orthogonal to ipj- 

If we denote by Tj the temperature of the thermal bath formed by ma- 
jorons and majorinos and by the one relative to the "standard model" 
particles, for < T^d we do have Tj = a(T^)T^ where 



model" particles and g^:s(T^) counts the effective relativistic degrees of free- 




(10) 




(11) 



takes into account the various annihilation thresholds for massive "standard 
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dom contributing to the entropy density. 

It is easy to show that majorinos can annihilate into a pair of majorons 
only through a p-wave and denoting the thermally averaged cross section 
for the process 'ipj'ipj J J by (crlfl) = gq x (Tj/m^j), we can calculate 
the temperature T* at which majorinos freeze out from the thermal bath of 
majorons by simply comparing the interaction rate 
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m^j /Tj 



ao X 



(12) 



with the expansion rate of the Universe. Here we do not make use of the 
Boltzmann equation for n^j since it would give a result different from ours 
only for a few percent. The freeze out temperature is then 
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where 



A 



a{T; 



1 1"^ 
\nA - -In (InA) 



1/2 



0.038 m^^Mpi 



(13) 



(14) 



9* [t; 

The today contribution of majorinos to the Q parameter is then given by |[I2| 

3/2 



2.4 X 10^*^ 



2tiT* 



X a 



/ T-i* \3/2 rpZ 

K-'-'y) -'7, 



tod 



'y,tod 



GeV" 



(15) 



9*s{T*] 

where T^^tod — 2.75 K is the today temperature of the relic photons. 

Taking h ~ 0.5 and ~ 0.6 as suggested by the latest C + HDM 
simulations and a(T*) ~ .5, we get ao — 10^^^ GeV~^. Since a detailed 



calculation of ctq with standard techniques [TTI] gives, for A ~ 0.5, ctq — 2 x 
10~^(m^^/V^^) |13|, we obtain that majorinos can form the CDM component 
of the Universe for a relativley small value of Vbl, Vbl — 10^ GeV. Such a 
small value of Vbl is quite likely to manifest itself in rare decay processes. 
However, in the model as presented, it would be rather difficult to detect the 
CDM component either through direct searches or through indirect detection 
of annihilation products of majorinos that annihilate in the Sun, in the Earth 
or in the galactic halo since the majorino is very weakly coupled to matter. 
It is however possible to add to the superpotential in Eq. (2) terms like 
ZHiHd + Z"^, which then allow the Majorino to have somewhat stronger 
interaction with matter, while still remaining a viable CDM. This will be 
discussed in more detail in a forthcoming publication |14 . 



Let us now turn to the nature of the neutrino spectrum. As it is well- 
known, the neutrino mass matrix in this model has the familiar see-saw form 
and leads to neutrino masses given by m,^- ~ (^^^^vt^if/^i)- The nature of 
the light neutrino spectrum depends crucially on the choice of h,y.^s, which 
are apriori arbitrary. If we choose the Yukawa couplings h^. ^ 10^^ (which 
is comparable to the Yukawa coupling of the electron and is therefore not 
unreasonable), then all three neutrino species have masses in the eV range 
and get to play the role of the hot dark matter together. On the other hand, 
it is perhaps more suggestive to have h^.vwK — rni-, in which case, one may 
get rrii,^ ^ 100 to 200 keV and mj^^ ~ 5 to 10 MeV. In our model, the heavier 
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neutrinos (i.e. z/^ and Ur) decay to the lighter neutrino plus a Majoron. Their 
decay is suppressed in the lowest order and arises in order ^ {jnyjVBLY- 



These decay rates have been calculated in detail in ref. [0, where it is found 
that for z/j Vj + J decay, the lifetime Tij is given by 



m,,.mf 



rr' ^ sin^2a 7 (16) 
^ 167r 

Here, a is a rotation angle in the Majoron coupling matrix and being 
unrelated to the observed neutrino mixing angles, its value can be of order 
one, although we will allow it to be anywhere from 0.1 to 1 in our discussion 
below. From the above equation we see that, for a 10 MeV tau neutrino, 
both the cosmological mass density fl^ as well the galaxy formation flS 



constraints are easily satisfied and the decay neutrinos , which get redshifted 
sufficiently to leave the universe matter-dominated before the epoch of galaxy 
formation, join the already present relic HDM neutrinos. Coming to the i/^ 
since its mass has to be less than 220 keV , it cannot satisfy either of these 
constraints. As a result, in our model z/^ must have a mass in the few eV 
range contributing to the HDM content. Thus, to summarize, the HDM can 
contain both the z/g and z/^ with masses in the 2 to 3 eV range each whereas 
the tau neutrino can either be very light and in the eV range or in the heavy 
range of several MeV's. In the latter case of course, it decays to both the 
electron and muon neutrinos. Also it is worth noting that, the existence of 
a Majorana z/g with mass in the 2-3 eV range will be testable in the ongoing 
PPou experiments. Failure to observe such a mass will imply that in our 
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picture, a few eV muon neutrino becomes the dominant HDM. 

In conclusion, we have given an example of an extension of the MSSM, 
where it is possible to implement the C + HDM scenario for large scale 
structure formation using a single new scale for physics beyond the standard 
model, i.e.the scale of U{1)b-l breaking scale needed anyway to understand 
the possible small neutrino masses. It is interesting that, the relative abun- 
dances of the CDM and HDM components are set by the same scale being in 
the TeV range, raising the hope that such theories may be testable in near 
future. 
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